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ABSTRACT 

There is accumulating evidence that at least a fraction of binary neutron star 
mergers result in rapidly spinning magnetars, with subrelativistic neutron-rich 
ejecta as massive as a small fraction of solar mass. The ejecta could be heated 
continuously by the Poynting flux emanated from the central magnetars. Such 
Poynting flux could become lepton-dominated so that a reverse shock develops. 

It was demonstrated that such a picture is capable of accounting for the optical 
transient PTFllagg (Wang & Dai 2013b). In this paper we investigate the X- 
ray and ultraviolet (UV) radiation as well as the optical and radio radiation 
studied by Wang & Dai (2013b). UV emission is particularly important because 
it has the right energy to ionize the hot ejecta at times t < 600 s. It is thought 
that the ejecta of binary neutron star mergers are a remarkably pure sample 
of r-process material, about which our understanding is still incomplete. In this 
paper we evaluate the possibility of observationally determining the bound-bound 
and bound-free opacities of the r-process material by timing the X-ray, UV, and 
optical radiation. It is found that these timings depend on the opacities weakly 
and therefore only loose constraints on the opacities can be obtained. 

Subject headings: radiation mechanisms: non-thermal — stars: neutron 


1. Introduction 

Compact binary mergers are of great astrophysical importance as the primary sources 
for the upcoming next generation ground-based gravitational wave (GW) detectors (see 
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Abadie et al. 2010; Bartos et al. 2013, for recent reviews). The detection of GW can be 
confirmed if accompanying electromagnetic (EM) signals are detected at the same time. 

It has been believed that the mergers of compact binaries, i.e., donble nentron stars 
(NSs) or a NS with a stellar-mass black hole, are the progenitors of short gamma-ray 
bnrsts (SGRBs; Paczyhski 1986; Eichler et al. 1989; Barthelmy et al. 2005; Fox et al. 2005; 
Gehrels et al. 2005; Rezzolla et al. 2011). The evidence for snch a scenario has recently 
gained strong snpport when Berger et al. (2013) and Tanvir et al. (2013) discovered an r- 
process kilonova associated with the SGRB 130603B^. Kilonova, powered radioactively by 
the nentron-rich material ejected dnring the coalescence of compact binaries, was hrst an¬ 
alytically predicted by Li & Paczyhski (1998) and then intensely stndied by many anthors 
(Knlkarni 2005; Rosswog 2005; Metzger et al. 2010b; Roberts et al. 2011; Metzger &: Berger 

2012) . Other EM signals inclnde radio afterglows (Nakar & Piran 2011; Metzger & Berger 
2012; Piran et al. 2013; Rosswog et al. 2013) and possible X-ray emission prodnced by the in¬ 
teraction of the NS magnetospheres dnring the inspiral and merger process (Palenznela et al. 

2013) . 

In a popnlar scenario, the hypermassive remnant of binary nentron star (BNS) merger 
collapses into a black hole on a timescale of ~ 200 ms (Faber & Rasio 2012). Models involv¬ 
ing a black hole snrronnded by a hyper accretion disc have been proposed (Popham et al. 
1999; Narayan et al. 2001; Kohri & Mineshige 2002; Xue et al. 2013; Kawanaka et al. 2013) 
as the central engine of GRBs. Several anthors (Dai et al. 2006; Zhang 2013), however, 
snggest that a highly-magnetized, rapidly spinning stable NS (magnetar) conld form dnr¬ 
ing the coalescence of BNSs. The evidence for magnetar formation following some SGRBs 
is accnmnlated dnring the past decade (Dai et al. 2006; Fan &: Xn 2006; Gao & Fan 2006; 
Rowlinson et al. 2010, 2013; Dai & Lin 2012; Wang & Dai 2013b; Gompertz et al. 2014). In 
addition, magnetar activity conld also be responsible for the statistical properties of X-ray 
flares from some GRBs (Wang & Dai 2013a). Recent simnlations in nnmerical relativity 
(e.g. Giacomazzo & Perna 2013) conhrmed that formation of a stable long-lived magnetar is 
a plansible scenario. However, in some cases the presence of snch stable long-lived magnetars 
has been rnled ont by radio observations, as in the cases of GRB 050724 and GRB 060505 
(Metzger & Bower 2014). 

Zhang (2013) recently proposed that if hypermassive millisecond magnetar resnlting 
from a coalescence of BNSs snrvives for snfficiently long time, a signihcant fraction of NS-NS 


^Jin et al. (2013) put forward an alternative scenario based on a two-component jet model. The non¬ 
detection of late radio emission of GRB 130603B (Fong et al. 2014) does not favor such an interpretation 
though. Despite of this fact, other possibilities (Fan et al. 2013; Takami et al. 2014) cannot be ruled out on 
current observational ground. 
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mergers would be discovered as bright X-ray transients associated with GW bursts without 
SGRB association. Gao et ah (2013) further studied the broadband EM signals of forward 
shock (FS) driven by the ejecta based on the energy injection scenario (Dai & Lu 1998a,b; 
Zhang & Meszaros 2001). Subsequently, Yu et ah (2013) studied the supernova-like EM 
signals of ejecta (dubbed merger-nova) powered by the Poynting flux from post-merger mag- 
netars by assuming the complete absorption of Poynting flux by ejecta. Wang & Dai (2013b), 
on the other hand, analytically studied the reverse shock (RS) emission powered by the cen¬ 
tral millisecond magnetars based on the fact that the Poynting flux is more likely to become 
lepton-dominated (e’''e“ pairs; Goroniti 1990; Michel 1994; Dai 2004; Yu & Dai 2007). Such 
a scenario successfully interprets observational data of the optical transient PTFllagg dis¬ 
covered by the Palomar Transient Factory group (Genko et ah 2013; Wang & Dai 2013b)^. 
More recently Metzger & Piro (2014) investigated the effects of pair-creation and annihila¬ 
tion on the optical and X-ray emission from the remnant of BNSs mergers powered by stable 
millisecond magnetars. 

One of the major concerns in the theory of BNS mergers is the uncertainty in the optical 
and nuclear properties of the ejected matter, which have started to attract attention only 
recently (Goriely et ah 2011; Korobkin et al. 2012; Barnes & Kasen 2013; Bauswein et al. 
2013; Grossman et al. 2014; Tanaka & Hotokezaka 2013). In particular, the opacity of the 
ejecta was usually set as k = 0.2cm^g“^ in the previous studies (e.g., Li & Paczyhski 1998; 
Rosswog 2005; Roberts et al. 2011; Yu et al. 2013), which was recently found to be k ~ 

10 cm^ g“^ (Barnes &: Kasen 2013; Kasen et al. 2013; Tanaka & Hotokezaka 2013; Grossman et al. 
2014). As a result of much higher values of opacity, much redder and dimmer kilo nova tran¬ 
sients are expected, with their peak shifted to infrared, just as in the case of event associated 
with SGRB 130603B (Berger et al. 2013; Tanvir et al. 2013). 

The ejecta from a neutron star merger (NSM) are thought to be a pure sample of r- 
process material. Apart from the above theoretical studies of the r-process material, it will 
be particularly valuable to observationally constrain the opacity of the r-process material. 

In this paper, we evaluate the possibility to constrain the bound-free and bound-bound 
opacities, which are strongly dependent on the nuclear composition of the ejecta. The 
physical picture of this scenario is illustrated in Figure 1. High opacity of the ejecta is 
caused by the presence of elements with half-hlled /-shells, such as lanthanides and actinides 
produced in the r-process (Kasen et al. 2013). As a result of high (bound-bound) opacity, 
optical emission from RS would be blocked during the early expansion of the ejecta. Early 
X-ray emission from the RS, on the other hand, would be blocked only by the electron 


^Wu et al. (2014) proposed an alternative interpretation for this optical transient within the framework 
of a magnetar-powered blast wave. 
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scattering opacity and unaffected by the bound-bound transitions. Consequently, we can 
determine the launch time of the Poynting flux more accurately, more dehnitively by timing 
X-ray emission. UV radiation from RS, however, is of great importance because it has 
the right energy to ionize the hot ejecta. Through a fortunate observation of X-ray, UV, 
and optical emission from RS powered by a millisecond magnetar, we can determine the 
launch time of the central magnetar (by X-ray observation), ionization break out by the UV 
radiation and the time when the ejecta become optically transparent. The later two times 
depend on the bound-free and bound-bound opacity, respectively. 

In Wang & Dai (2013b), we demonstrated how our model consistently accounts for the 
observed optical and radio properties of the PTFllagg transient. We considered the case in 
which the spin-down time of the magnetar Tgd exceeds the deceleration time of the blast wave 
Tdec (e.g.. Case I of Gao et ah 2013). In Section 2, we further study the other two cases. In 
our calculations, we adopt an improved version of the blast wave dynamics, which is different 
from the one used in Gao et ah (2013) and Wang &: Dai (2013b) in that it accounts for the 
different regions in the blast wave. We compare different prescriptions for the dynamics in 
Section 2. In the same section we additionally investigate the effect of the optical obscuration 
by the r-process material, neglected in Wang & Dai (2013b). We study the ionization effect 
of the UV radiation on the ejecta in Section 3. Finally, we discuss how to observationally 
constrain the bound-bound and bound-free opacities in Section 4. A summary is given in 
Section 5. 


2. Leptonized reverse shock emission 

Numerical simulations suggested that the ejecta from an NSM have a typical velocity 
V = 0.1-0.3c and mass Mej = Mq (Rezzolla et al. 2010; Hotokezaka et ah 2013; 

Rosswog et al. 2013). The resulting central compact object, under the assumption that 
it can avoid collapse and persist for a long time, would dissipate its rotational energy by 
launching a luminous Poynting flux. Without the impact of the Poynting flux from the 
nascent millisecond magnetar, the matter ejected during the merger would have spatially 
extended morphology, which would soon have started to expand homologously, despite the 
effects of radioactive heating (Rosswog et al. 2014). With the onset of Poynting flux, the 
ejecta can be rapidly compressed into a thin shell by the overpressure caused by the energy 
injection. The overpressure will lead to a shock driving into the ejecta. The postshock 
pressure is p 2 = (7 + 1) for strong shock (Landau & Lifshitz 1987), where Vi is 

the unshocked fluid velocity relative to the contact discontinuity, pi its density. With pi = 
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Mej P2 = Lsdt [3 X Itt {v + and the width of the initial ejecta A = vit, 

we hnd the crossing time 


_ A _ 

tnrnss 


67 A^MejV^^ 


Vi V7 + I L, 


sd 


= O^OSsAf 


( 1 ) 


and fi at this time 


wi, cross = 3.8 X < n ^ 0.2c, 


( 2 ) 


where Lsd is the spin-down luminosity of the magnetar. Here the usual convention Q = 
is adopted. This crossing time is negligible compared to the activity duration of the central 
magnetars, justifying the model we are considering. 

As mentioned in Section 1, the Poynting flux from millisecond magnetars is more likely 
to become lepton-dominated. In developing a model to account for the transient source 
PTFllagg, Wang & Dai (2013b) adopted a dynamics 


Lq min {t, Tsd) = (7 - 7 ej,o) ^ejC^ 2 ( 7 ^ - l) 


( 3 ) 


which is different by a factor of 2 on the second term from that used by Gao et ah (2013). In 
the above equation Lq = ^Lsd with ^ the fraction of the Poynting flux catched by the ejecta, 
7 the Lorentz factor of the ejecta with initial Lorentz factor 7ej,o, the mass swept up by 
the shock. 


In this paper we would like to adopt an alternate dynamics similar to that used by 
Yu et ah (2013) and to evaluate any differences between these two kinds of dynamics. In 
this scenario, the ejecta gain energy because of the work done by RS and likewise the FS 
gains energy by the work done by the ejecta. As a result, the ejecta are sandwiched between 
RS and FS. The treatment of RS is otherwise similar to that in Wang & Dai (2013b). 

The total non-rest energy of the ejecta and shocked media (including FS and RS) can 
be expressed as 

Ek = -fE'^ + (7 - 7ej,o) MejC^ + 7^ej,int + “ l) , (4) 

where E'^ and E'^- are the respective energies of the reverse-shocked wind (region 3) and 
ejecta in the comoving frame. For the dehnition of regions 1-4, see Figure 1 (see also Dai 
2004). Here we neglect the rest energy of region 3 because it is lepton-dominated and 


^Here we ignore the initial size of the ejecta. 
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therefore very hot. In a way similar to the generic dynamic model for GRB afterglow 
(Huang et ah 1999), the dynamics can be derived as (cf. Yu et ah 2013) 


d'y 

dt 


^-hsd + L 


rd 


^ej,e 




dE. 


ej ,int 


dt' 


(^2 _ ^2 


Eg + MejC^ + £'ej,int + 


( 5 ) 


where V is the Doppler factor, the radioactive luminosity of the ejecta in the observer 
frame, Lej,e the energy loss rate in ejecta due to thermal radiation. The evolution of energies 
in region 3 and ejecta can be expressed as 


dt' 


dE'^ 


ejjint 


dt' 


dVU, 


T f 3,giic ( r '\ 

dV^ dV- 

^*^3,enc ^’^ej,enc , rf t' /^\ 

where t' is the time measured in the comoving frame. Lb g the thermal radiation of ejecta. 
Here we have to discern different volumes involved in this problem. See Figure 1 for reference. 
Vjj'g^g and lY gjjg are the comoving volumes enclosed by region 3 and ejecta, respectively. 
and 1/j, on the other hand, are the volumes actually occupied by these two regions. In other 
words, the volumes with the subscript enc denote the spherical volumes, the volumes without 
such subscript denote the shell volumes. The pressures in region 3 and ejecta in the comoving 
frame are = e^/S and Pej = -^ej,int/3Kp respectively. The width of the sandwiched ejecta 
can be determined by assuming pressure balance ps = pey Because the ejecta are sandwiched 
in a thin shell, i.e., the width of ejecta (in the comoving frame) Ab <C r, so that its volume 
0, Equation (7) can be approximately written as 


V- 

ej 


dt' 


— T' — T' 

— -^rd -^ej,e- 


The comoving volume can be found as (Yu et al. 2013) 

dV 

= Airr^^fdc, 


with 


dt' 

dr 

dt 


fdc 

rws' 


( 8 ) 


( 9 ) 


( 10 ) 


The radioactive luminosity in the comoving frame = Lrd/F^^ is determined according to 
Equation (4) of Korobkin et al. (2012). The thermal energy Lb ^ emitted by the sandwiched 
ejecta in the comoving frame is still expressed by Equation (7) in Yu et al. (2013), while the 
luminosity light curves at a particular observational frequency are modihed as 


vL,, = 


Ab 


{hv/Vf 


max (r, 1) r exp {hu/VkT') — 1 ’ 


( 11 ) 
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i.e., there is an extra factor A'j/r compared with the solid spherical geometry. 

Inspection of Equation (3) indicates that initially (7 — 7ej,o) ^ 2 ( 7 ^ — 1) 

the ejecta will be accelerated linearly in time until t = min (Tgd, Thee), where the deceleration 
timescale T^ec is determined by the condition (7 — 7ej,o) AfejC^ = 2 (7^ — 1 ) (see also 

Gao et al. 2013). By setting the spin-down timescale Tsd ~ T^ec, we will arrive at a critical 
ejecta mass (Gao et al. 2013) 

Mej,, ~ (12) 

where I is the moment of inertia of the neutron star. It is this critical mass that dehnes the 
three cases considered by Gao et al. (2013), i.e., Gase I for Tgd > Tdec, Gase II for T^d = Tdec, 
and Gase III for Tsd < Tdec- Here we show the corresponding results including RS in Figures 
2-4, where the spectral index of electrons in FS is set as = 2.3, and the luminosity 
distance of the source is set as = 10^^ cm. The optical opacity is set as k = lOcm^g”^. 
For more discussions about the dependence of opacity on wavelength, see Section 3. To ease 
comparison between different choices of dynamics, we selectively show in Figures 2a and 2b 
as dashed lines the results with the same parameters as the solid lines but with dynamics 
expressed by Equation (3). 

In the above calculations we set the Lorentz factor of the unshocked electrons/positrons 
(region 4) as 74 = 10^, as determined in the literature (Atoyan 1999; Dai 2004; Wang & Dai 
2013b). The spectral index of e+e“ in RS is set as = 2.2 (see, e.g., Wang & Dai 2013b). 
There is an additional complication concerning RS that should be mentioned. After the 
shut-off of the central magnetar at Tsd, the reverse shock crosses region 4 so that there are 
no more e'''e“ to be shocked. Gonsequently, the electrons already cooled to a Lorentz factor 
7c cannot be shocked again to a higher Lorentz factor. As a result, in the analytical and 
numerical calculations, if we hnd the cooling Lorentz factor increases in a time period after 
Tsd, which is the case when t > T ]^2 for Gase I, we will set its value the same as before. 

We hnd that the characteristic frequencies of synchrotron radiation of FS are quite 
similar to that calculated by Gao et al. (2013) so that we do not show them in the hg- 
ures. Because the RS emission will be absorbed by the ejecta, we show the RS emission 
with/without absorption in Panels (d) and (c) respectively. It is in principle possible that a 
FS and a RS develop in the ejecta upon the interaction with region 3 and region 2 respec¬ 
tively. In practice, however, since the ejecta are very thin, the FS and RS in it can last only 
for a transient while. Gonsequently we do not consider them here. 

A hrst glimpse of Figure 2a shows that the evolution of Lorentz factor with different 
choices of dynamics are very similar. However, there are also differences that result in 
appreciable modiheation to the evolution of characteristic frequencies (comparing the solid 


lines and dashed lines in Figure 2b) and light curves, which are not shown in Figure 2. The 
reason is that the ejecta absorb some amount of energy, resulting in a heavier ejecta and 
therefore a delay of the deceleration time T^ec by a factor of 2.1. This further results in a 
drop of Um (see Figure 2b), which is otherwise identical to the situation where the ejecta do 
not absorb energy (comparing the solid lines and dashed lines in Figure 2b). The delay of 
Tdec also strengthens the FS since more ambient media are swept up at T^ec- Figure 4a shows 
the coasting of Lorentz factor, but not so clear-cut as in Figure 4a of Gao et ah (2013). 

Although the energy absorption by ejecta has an appreciable effect on the dynamics and 
therefore the light curves, it is not so signihcant as for the case of merger-nova discussed by 
Yu et al. (2013). We therefore conclude that a simple model based on Equation (3) and used 
in Wang & Dai (2013b), properly reproduces qualitative features of the more sophisticated 
one employed here, and therefore the use of this model in Wang & Dai (2013b) is justihed. 


For completeness and also for the ease of future quantitative analysis, we present the 
analytical results of RS based on the dynamics (3) for Case II and Case III below and their 
temporal scaling indices in Table 1 (the corresponding results for Case I can be found in 
Wang & Dai 2013b). For Case II, the various time scales and the peak Lorentz factor are 


Tm 

Tct 

Tac 

Tdec 

Tm 

T 

am 2 


= 2.1 X 10 '^daysLo, 49 ^ej ,-4 
= 2.8 X 10-3 daysL- 49 '^e/- 44 (-i 

^ Tmc = 5.0 X lO-^daysL-^/^Mj/I^eJ/J.ieyM// 
^ Tsd 

= 2.4 X 10^ daysLl%T^J^ln-T^ 

= 9.4X 


X 


(p-i)r 


3 p + 22 


12 


3 p +2 
12 


P 


p 


25 


7sd — 28Lo,49^sd,3Afej^^4 -|- 1, 


(13) 

(14) 
(16) 
(16) 
(17) 


(18) 

(19) 


where Tam 2 is the time when crosses at the second time, which is very sensitive to many 
parameters because these two frequencies almost have the same temporal evolution indices 
during the time period Tm < t < Tam 2 (Table 1, see also Figure 3b). The characteristic 
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frequencies and the observed peak flux are 


^a,sd 

1.7 X 

V /\/f'^(P+2)/(P+^)^:2(p—l)/(p+4) (p+2)/2(p+ 4) 2(p—2)/(p+4) 
^^‘^ej,-4 74,4 

(20) 

^m,sd — 


(21) 

^c,sd 

1.9 X 

(22) 

-fl^,max,sd 

3.5 X 10^mJyLQ4g . 

(23) 


For Case III the corresponding values are 


and 


^N1 

= 2.1 X 10 ^daysLg jgMej _3 

(24) 

Fdec 

= 0.9 daysL„-I'=’T.7f 

(25) 

Tm 

= 89.1daysLJ2>r75ir‘''“ 

(26) 

7sd 

= 5.6Lo,49Tsd,3Mej^^3 -|- 1, 

(27) 


^a,sd 

= 6.9 X 10^3 




, .2(p+5)/(p+4) (p+2)/2(p+4) -2/(p+4) 

^B-1 /4,4 


( 28 ) 

^m,sd 

= 2.3 X 10 “ 

t )^ 

( 29 ) 

^c,sd 

= 1.8 X 


( 30 ) 

■^^',max,sd 

= 4.6 X 10 ^ JyLg 4g _3e^^_474 4^27 . 


( 31 ) 


3. Ionization Break Out 
3.1. Opacity Estimate 

Our calculations above are based on the results of Kasen et al. (2013), who studied 
the opacity at optical/infrared wavelengths, in particular the bound-bound opacity. For 
UV radiation, we are interested in the radiation at early times, i.e. t < 600 s, because the 
ejecta will become transparent at later times. During such early times, the ejecta are hot 
enough to thermally ionize the lanthanides, e.g. Ce, to Ce^+. For simplicity, in the following 
calculations we assume that Ce is completely thermally ionized as Ce^+. As estimated below. 
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the opacity at UV wavelengths is dominated by the bonnd-free (photoionization) transitions. 
We will therefore calcnlate the photoionization of Ce^+ by the UV radiation from RS. The 
forth ionization potential of Ce is X 4 = 36.72eV (Cox 2001). We hrst estimate the different 
opacities as follows. 

The wavelength independent electron scattering opacity is (e.g., Kasen et ah 2013) 

where A is the mean atomic weight of the ions, x the mean ionization fraction, nip the proton 
mass. After nltraviolet ionization, Ce is ionized to a level x = 4. With A ^ 140 we hnd 
Kes ~ 0.01 cm^ g“^, which is a factor of ~ 20 smaller than the nsnal valne Kes = 0.2 cm^ g“^. 

The free-free opacity can be fonnd as (e.g., Rybicki & Lightman 1979; Kasen et ah 2013) 

= O.lS^p-gTs-'/'A^ 5 (1 - cm^ g-^, (33) 

where the typical density of the ejecta at early times is p > 10“®gcm“^ and temperatnre 
T >2 X 10^ K ~ Tion = Xs/fc, where xs = 20.198 eV is the third ionization potential of Ce, 
and Tion the corresponding temperatnre. For the valnes x and A mentioned above, we see 
that the free-free opacity at UV wavelengths is Kg ps 5 x 10“®cm^g“^, which is completely 
negligible. The free-free opacity at X-ray band is even smaller because of the dependence of 
opacity on wavelength A. 

The bound-bound opacities of the lanthanides are a function of both temperature and 
wavelength (Kasen et ah 2013). On the one hand, the line expansion opacity of cerium (Ce, 
Z = 58, /-shell), for example, decreases to the red and sharply drops to zero for A < 1000 A 
(Figure 7 in Kasen et ah 2013). That is, the bound-bound opacity in the UV and X-ray 
bands is negligible at the temperature we are interested in, i.e. T > Tjon- On the other hand, 
the Planck mean expansion opacity of neodymium (Nd, Z = 60, /-block, similar to Ce in 
shell structure) increases sharply with temperature when T < 4000 K and then decreases 
with temperature (Figure 6 in Kasen et ah 2013). The reason is that, with the increase 
of temperature, more excited levels are populated. But if the gas becomes hot enough to 
ionize, the originally populated state leaves blank and the opacity declines. The bound- 
bound opacity declines sharply when the temperature T > 15000 K, i.e., Kbb ^ 10“^cm^g“^ 
(Figure 9 of Kasen et ah 2013). For the temperature we are interested for UV radiation, i.e. 
T > Tion, becomes negligibly small. 

The bound-free opacity at the threshold energy is given by (e.g., Kasen et al. 2013) 

(To 

“ im, Z (T) 


( 34 ) 
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where Z (T) is the partition function and Ai? the excitation energy. The Boltzmann factor 
^-AE/kT account of the fact that an atom has to be thermally excited so that a photon 
has enough energy to ionize the excited atom. The wavelength-dependent parameter ao is not 
accurately constrained but can be approximated as the hydrogenic value (Tq ~ 6 x 10“^^ cm^ at 
optical/UV wavelengths. For UV radiation, AE ^ 0, unless the temperature T ^ Tion- The 
largest uncertainty for the calculation of Kbf comes from the partition function Z (T). Irwin 
(1981) tabulated Z (T) for Ce'*' and Ce^'*' in the temperature range 1000 K-16000 K. The 
partition function increases slowly with temperature. It is also expected that the partition 
function of Ce^'*' is smaller than Ce^"*" at high temperature. Accordingly, we estimate Z (T) ~ 
1000 for T > Tion- Taking the typical values, we hud ^ cm^ g“^ for T > Tjon- 

We therefore conclude that for UV photons, bound-free opacity is dominant. For X-ray, 
AE = 0, but (To takes a value several orders of magnitude smaller than it does at the 
optical wavelengths (Hakken & Wolf 1987). Consequently, the bound-free opacity for X-ray 
is negligible and the opacity at X-ray is dominated by Kes- 

Although the above reasoning is quite robust, but owing to the complex nature of 
opacities, it is worth of more scrutiny here. Figure 7 in Kasen et ah (2013) shows that 
the opacity of osmium (Os, Z = 76, ci-shell) is very high at UV wavelengths. But this 
high opacity is only for temperature ~ 5000 K. Figure 6 in Kasen et ah (2013) shows that 
the opacity of Fe (also d-shell in shell structure, therefore similar to Os) keeps constant 
in the temperature range 5000-14000 K. But with a temperature high enough to ionize 
Os, the opacity will decline, as explained above. Consequently we can tentatively infer 
that the bound-bound opacity of Os at temperature Tion is smaller than the bound-free 
opacity evaluated in the previous paragraph. Another uncertainty comes from actinides 
(90 < Z < 100). While the actinide series is generally of very low abundance, but because 
of the fact that line expansion opacities are expected to shift to shorter wavelengths for 
heavy elements such as actinides, the actinide series may make signihcant contribution to 
the opacity that exceeds the bound-free opacity estimated above. If the bound-bound opacity 
contributed by Os and/or actinides is in excess of the bound-free opacity, we will not observe 
an UV ionization break-out. 


3.2. Ionization by nltraviolet radiation from RS 

Before elaborating the calculation of ionization, we need evaluate the recombination 
factor, which can be calculated according to the bound-free opacity as (Rybicki & Lightman 
1979) 
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Or 


= 47r f 




\2TrkT / 


>/-Ei(^). (35) 

/ mtc^ QeQj^ Z (1 ) \k1 / 


/3 

h 


where (^e, 9+, On are the statistical weight factors of electron, Ce'^’*' and Ce^’*', respectively, 
rrie the electron mass, and Ei (x) is dehned as 


r-oo g-t 


Ei (x) = / — dt. 

.1 t 


(36) 


For Ce"^+ we set x = Xi- 


The width, Ah^, of the ionized zone (Ce'^'*') of the ejecta shell in the comoving frame 
evolves as (cf. Eqnation (9-8) of Harwit 2006) 


= (47rr2^n3+) ^ ^ - Ak ^Uearec, 


dt’ 


where 




n3+ = 


(37) 


(38) 


dTrr^^AgjArrtp 

is the nnmber density of Ce^’*' before ionization by UV radiation and dNi/dt' the injection 
rate of the ionizing photons from RS. The electron density rif. is rif. = xn 4 +, where 77 . 4 + 
( 77 , 4 + = 77 , 3 +) is the nnmber density of Ce'^’*'. By setting 774 + = 773 +, we neglect the thermal 
expansion of the photoionized zone due to the increase of pressure relative to the Ce^+ zone. 

The UV photons from RS are scattered off free electrons before ionizing Ce^+ zone. The 
optical depth for the ionizing UV photons is 


Tuv = 71 mpn 3 +K«'^ej,i' 


(39) 


The injection rate dNi/dt' is therefore evaluated as 


m 

dt' 


AttD 


p-njv 

2 


hu' 




(40) 


where F/, = V is the energy flux from RS in comoving frame, hv' = Xa is the energy of 
ionizing photons. The thermal broadening is given by 


Au' = 


3kT 

ArUpC^ 


(41) 


The thickness of the transition from Ce^"*" zone to Ce^’*' zone is of the order of mean free 
ionizing path <5 = ( 773 _|_(Tbf)~^ — 1.5 x 10® cnmij^’; ^ 4 , which is more than one order of magnitude 
smaller than the width of the ejecta Ah > 10 ^ cm. 
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The calculated UV light curves and ionization fraction (Agj j/Agj) are depicted in Figure 
5. The observed UV flux is attenuated by a factor where Tbf = hm (Ah — Ah j) Ampn^+. 
Before the ejecta are completely photoionized, Tbf > 100 so that effectively the UV flux can 
be observed only when the ejecta are completely photoionized, i.e. ionization break out. 
This is why we see an abrupt increase of UV flux when the ionization fraction equals 1. As 
the ejecta expand, the ejecta temperature drops rapidly so that when T = Tion the thermally 
ionized Ce^’*' begins to recombine. We see from Figure 5 that it is not long when the ejecta 
become transparent after the temperature drops below Tion- As a result, our treatment of 
the ultraviolet flux is not affected very much by the recombination. 


4. Discussion 

T-process material is of great importance astrophysically, the studies of which are un¬ 
fortunately very difficult in laboratory. It is therefore particularly valuable if we can study 
their properties by observing radiation, obscured by a pure r-process material, at different 
wavelengths. To this aim, we propose in this paper to observe the RS emission obscured by 
the ejecta, which is believed to be a pure r-process material. 

Figures 2-4 show that the ejecta shell is subrelativistic when it becomes transparent 
(at time tx.thin) to X-ray emission from RS. During this phase the radius of the ejecta shell, 
suitable for Case I, II, and III, is 

r = 2.1 X 10^=" t < Tni (42) 

from which fx,thin can be estimated as 

tx,thin = 7.1 X 10^ (43) 

upon substituting Equation (32). Inspection of Figures 2-4 indicates that tx,thin is approx¬ 
imately the time when the X-ray flux peeks from the left edge of the light curve. This 
behavior is also true for optical flux. Note that this estimate for fx,thin is accurate enough in 
comparison with the numerically determined value. 

The estimate for the time, topt,thin, when the ejecta become optically thin for visible 
light can be done in a similar way as long as we are aware of the fact that at this time the 
ejecta shell is in relativistic regime so that the radius should be given by 

r = 1.3 X 10^^ cniTg 47 Mgyi 4 t 3 , Tati <t < min (T^d Tdec) (44) 

V.thin = 2.2 X 


from which we get 


(46) 
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Unfortunately, topt,thin depends on Kbb very weakly so that only a loose constraint on Kbb 
can be obtained. If, on the other hand, the ejecta are massive enough and/or the spin- 
down luminosity of the central magnetar is moderate so that the ejecta shell is still in the 
subrelativistic regime when it becomes transparent for visible light, we have topt.thin ot 
and Kbb can be constrained more compactly. 

To evaluate the dependence of tuv.nO) the break-out time of UV radiation, on Kbf, we 
vary Kbf by a factor /, i.e., = /fi^bf- From Equations (34) and (35) we see that the factor 

/ contains the uncertainties in (Tq and Z (T). The uncertainty in (Tq is a few, so is Z (T). 
So we expect the uncertainty of / is ~ 10. The solid lines from left to right in the insets 
of Figure 5 show the UV light curves with / = 0.1, 1, and 10, respectively. This time lies 
between fx,thin and topt,thin, as it should be. 

This model contains many free parameters to be constrained by htting light curves. 
Wang & Dai (2013b) demonstrated that by simultaneously htting the optical and radio light 
curves, one can accurately constrain the spin-down luminosity Lgd, the spin-down time scale 
Tsd, the ejecta mass Mej, the density of ambient medium (region 1) n, the Lorentz factor of 
leptonized wind (region 4) 74 , the electron power-law index of reverse shocked wind (region 3) 

the magnetic energy fraction and the redshift 2 of the source (Table 2 of Wang & Dai 
2013b). 

Wang &: Dai (2013b) also determined the launch time of the Poynting hux of PTFllagg. 
This is possible because the launch time affects the temporal decline index of the early optical 
light curve. To accurately determine the launch time of the Poynting hux by timing X-ray 
hux. Figures 2-4 indicate that we have to take observation in X-ray band following the merger 
no later than ~ 100-700 s. A similar observational strategy has to be employed to determine 
the ionization break-out of ultraviolet light (see Figure 5). Future gravitational observation 
of NSMs may also help determine the launch time of Poynting hux because this time roughly 
coincides with the merger time. Figure 5 shows that the ionization break-out time depends 
on Kbf weakly. Thus by timing X-ray, UV, and visual light we can obtain weak constraints 
on Kbb and Kbf- 


5. Conclusions 

In this paper, we have analyzed electromagnetic counterparts in various bands to the 
binary neutron star merger, and arrived at the following conclusions: 

(i) By assuming that the Poynting hux of the central magnetar becomes lepton-dominated 
so that RS develops, broad-band EM signals, i.e. radio, optical, UV, and X-ray emis- 
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sion, can be produced. The NSM ejecta can be accelerated to relativistic speed so that 
a forward shock can form. We therefore expect to observe EM signals hrst from RS 
and subsequently from FS. 

(ii) In the study of the optical and radio radiation of PTFllagg, the simple dynamics, 
i.e., Equation (3), was adopted by Wang & Dai (2013b). In this paper, the dynamics 
is determined by pressure balance and it is found that Equation (3) is an acceptable 
approximation. 

(iii) In studying PTFllagg, Wang & Dai (2013b) ignored the absorption of optical radia¬ 
tion by the ejecta. We hnd in this paper, by including the absorption effect, that by 
the time the hrst optical datum was taken, the ejecta become optically thin so that 
the treatment of Wang & Dai (2013b) is justihed. Wang &: Dai (2013b) also ignored 
the emission from the ejecta itself, which is verihed in this paper (Figure 2) that the 
emission of the ejecta itself is negligible in all relevant bands. 

(iv) Scattering off free electrons will obscure the early X-ray emission from RS. Optical 
emission from RS will be suppressed by bound-bound opacity until the ejecta become 
transparent at later times. UV radiation from RS, on the other hand, has the right 
energy to ionize the hot ejecta during its early expansion and ionization breakout 
should therefore be observed. 

(v) Given the weak dependence of arrival times of EM signals at different wavelengths on 
different opacities, we conclude that the timing of X-ray, UV, and optical light is not 
an accurate method to observationally constrain the opacities of r-process material. 

(vi) Recent papers by Gao et ah (2013) and Wang & Dai (2013b) make an approximation 
that the merger ejecta is rapidly compressed into a thin expanding shell, similar to 
the one which forms at GRB afterglows. In this paper we estimate the time it takes 
for the flux to compress the ejecta to a thin shell and hnd that the time is negligible 
compared with the activity duration of the central magnetar. As a result, the models 
proposed by Gao et al. (2013) and Wang & Dai (2013b) can be applied safely. 


We are grateful to the referee for insightful comments and constructive suggestions, 
which signihcantly improve the presentation of this paper. This work is supported by the Na¬ 
tional Basic Research Program (“973” Program) of Ghina under Grant No. 2014GB845800. 
L.J.W. and Z.G.D are also supported by the National Natural Science Foundation of Ghina 
(grant No. 11033002), and Y.W.Y. by grant No. 11473008. 
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Fig. 1.— Schematic illustration of the physical picture described in this paper. The merger 
ejecta are catched up by the relativistic wind, driving the surrounding medium to form 
forward shock. The magnetar wind, on the other hand, is heated by the reverse shock. The 
ultraviolet radiation from the reverse shock ionizes the thermally ionized Ce^’*' to Ce'^’*'. 6 
is the transition width between Ce'^’*' zone and Ce^’*' zone. When the ejecta temperature 
drops below Tjon, viz. the temperature to thermally ionize Ce to Ce^’*', Ce^+ begins to 
recombine. Consequently, the ionization stage can only last for a limited time. The Arabic 
numbers indicate different regions: region 1 is the unshocked medium, region 2 is the shocked 
medium, region 3 is the shocked wind, region 4 is the unshocked cold wind. 
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Fig. 2.— Calculation results for Lsd,i = 10"^^ergs“^, Tgd = lO^s and Mgj = (a) The 

evolution of Lorentz factor, (b) The characteristic frequencies of RS emission. The three 
dashed lines mark the X-ray, optical (R) and radio (8 GHz) bands, respectively, (c) Light 
curves of RS without absorption by ejecta, (d) Light curves of RS with absorption by ejecta 
considered (thick) and of FS (thin). The two vertical dotted lines in panel (d) mark the times 
when the ejecta become transparent for X-ray and visual light respectively. The emission by 
the heated ejecta is negligible in this case, so we do not show them. Dashed lines in panels 
(a) and (b) are the results based on the dynamics expressed by Equation (3). In panel (a) 
Tni and Tn 2 are the times when 7 — 1 = 1 , i.e., the transition time between relativistic 
motion and Newtonian dynamics, Tct is the transition time for cooling Lorentz factor, see 
Wang & Dai (2013b) for more explanation. In the calculations we set initial ejecta velocity 

/^ej,0 = 0-2. 
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Fig. 3.— Calculation results for = 10'^®ergs = 10^s, Mej = 10 ^Mq. The 

meanings of the panels are the same as in Figure 2. Ejecta emission is negligible. 
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Fig. 4.— Calculation results for Lsd,i = lO'^^ergs”^, Tgd = 10^ s, Mej = The 

meanings of the panels are the same as in Figure 2. In panel (d) the optical emission by the 
ejecta is shown as thin line. Ejecta emission in other bands is negligible. This is the only 
case discussed in this paper that the ejecta emission becomes appreciable in optical band. 
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(a) Mej = 10 ^M0,Lsd,i = 10'^'^ergs ^ 




Fig. 5.— Ultraviolet radiation from RS. The solid lines correspond to the observed ultraviolet 
flux leaked from the ionizing ultraviolet flux from RS (dash-dotted lines). The dashed lines 
are the ionization fraction of the ejecta. The vertical dotted lines indicate the times when 
T /Tion = 1 and r = 1. The solid lines from left to right in the insets are the UV light curves 
by taking the bound-free opacity uncertainty factor / = 0.1, 1, and 10, respectively. 
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Table 1: Analytical Temporal Scaling Indices of Various Parameters of the RS. The Analytical 
Results for Case I were Presented in Wang & Dai (2013b). 


t < T] 


N1 


T, 


Ai < t < Tct 
Tct < t < T, 


ac 


sd 

V 2 


_ 7 - 1 ^ 

Case II: Lq ~ 10^^ erg ~ IQ-^Mq 

1 1 3p+14 33 

2 2(p+4) 2 2 

1 3 -5 -3 

p+4 

1 3 _5 9 

p+4 

1 3 -5 9 

p+4 

_3 1 9p+46 __9_ _U 

8 4 16(p+4) 16 16 

_6 2 3p+14 _3 _3 

5 5 5(p+4) 5 5 

_6 2 _]^ _3 _3 

5 5 25 5 5 

Lq ~ lO^^ergs”^, Mej ~ 10“^Mq 

1 3 3p+i4 33 31” 

^ 2 2(p+4) 2 2 

1 3 _32^ _5 _3 

P+4 

n 1 3p+14 _3 _3 

^ ^ 2(p+4) 2 2 

_3 1 9p+46 __9_ __9_ 

8 4 16(p+4) 16 16 

_6 2 3p+14 _3 _3 

5 5 5(p+4)5_5 


Tac<t< Z 
Tsd < t < Tj 
Tn2 < t < Tam2 
Tam2 7 t 

Case III: 


t < T/vi 
Tati < t < Tsd 
Tsd < t < Tdec 
Tdec < t < T/V2 

Tn 2 < t 


2 

-2 

-2 

-2 


16 

3 

'5 

3 

5 


1 

2 


-2 


3 

2 

_9_ 

16 

3 

'5 
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